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ABSTRACT 

A study of the fluctuating velocity field in the 
developing region of a pulsating pipe flow is presented. 

Theoretically, three separate determinations of radial 
and axial perturbation velocity profiles were calculated 
involving varying degrees of analytical simplification. 

Experimentally, pulsating disturbances were created 
by a sinusoidal pressure wave exciter installed at the out- 
let of a 1.5 inch circular wind tunnel facility. Cross- 
wire anemometer measurements of the pulsating velocity 
components compared favorably with the parallel theoretical 
results. Likewise, experimental determinations of the 
phase angles between velocity and pressure signals compared 
favorably with analytical values. 
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NOMENCLATURE 

Constants, coefficients of the polynomial 
approximation of the velocity profile 

Inside diameter of the pipe 

Pressure 

Radial. coordinate 
Radius of the pipe 
Reynolds number 
Time 

Mean velocity in the pipe 
Maximum velocity in the pipe 
Axial component of velocity 
Radial component of velocity 
Axial coordinate 
Kinematic viscosity 
Density 

Angle of phase shift 
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I . 



INTRODUCTION 



The development of flow in the inlet region of a circular 
cylindrical tube has been a subject of interest since the 
time of Hagen (1839) and Poiseuille (1841) . Osborne Reynolds 
(1883) conducted extensive investigations into the stability 
of pipe flow, and this area has become a classic problem 
with a sizable amount of literature published. However, the 
study of fluid flow through a rigid tube under the influence 
of a periodic pressure difference is more limited. 

The theoretical study of pulsating, laminar pipe flow 
was primarily given by Th. Sexl [1] and S. Uchida [2]. 
Experimental data are limited to the works of Atabek, Chang, 
and Fingerson [3] and the results of Florio and Meuller [4]. 
With the introduction of the modern computer, several attempts 
have been made to solve the problem by numerical methods . 

The most recent work in this regard is that of Lew and Fung 
[5], In the theoretical analysis the main problem has been 
the non-linear convective acceleration term in the Navier- 
Stokes equation. Uchida [2] was able to apply linear theory 
for the fully developed laminar pulsatinq flow. Lew and 
Fung [5] found this linearization to be valid in the entry 
region for high Reynolds numbers. However, this assumption 
leads to significant errors for small Reynolds numbers (the 
area of most interest) . 
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The present investigation was undertaken to study the 
fluctuating velocity field which arises in the developing 
region of a rigid walled conduit due to an imposed sinusoidal 
pressure gradient. Additionally, results were sought con- 
cerning any phase differences between the velocity and 
pressure signals with a view toward later stability 
determinations . 
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II. THEORETICAL ANALYSIS 



A. SOLUTION # 1 OF PULSATING FLOW EQUATIONS 

For oscillatory flow in a circular pipe, the Navier- 
Stokes equations must be solved for u and v, the fluctuating 
velocity components, as functions of three variables (r,x,t.) 
The velocities, u and v, are independent of <J> , the azimuthal 
coordinate. The applicable equations from Schlichting [6] 
are : 



The continuity equation is: 

+ Z = 0 ■ ( 3 ) 

3x 8r r 

It is assumed that u=u+u, p=p+p, and v = v where 
u and p are the steady components and u, v, and p are the 
pulsating components. 

After linearization, two additional assumptions are made 
namely, (1) that p is applied across the full cross-section 
of the tube, i.e., p is independent of r , and (2) that the 
gradient of v in the x direction is mild compared with the 
gradient in the r direction. 





( 1 ) 
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The resulting equations are: 
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Substitution of (6) 



into (5) yields 
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The left side of equation (7) is recognized as that treated 
by Uchida [2] and the right side of (7) is known if v and u 
are known. 

By assuming v = R (r) X (x) e"*" ^ t+< ^ , equation (4) may be 
solved by separation of variables. This technique yields 




for which the X solution is 
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Therefore the complete solution for v is 



v = C.a r e"“ x e i(Nt ' H,) 
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and -co is the constant introduced in the separation of 
variable technique. Now that a solution for v is known, 
equation (7) may be solved. By substituting equation (11) 
into (7) and further assuming 



u = b + b n r + b 0 r' 
o 1 2 



(14) 
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equation (7) becomes a zero order non-homogeneous Bessel 
equation of the form 

-int 
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(16) 



The relationship in equation (14) allows determinations in 
the developing region to be made. The relationships in 
equation (15) concerns an oscillating axial pressure gradient 
and a resultant oscillatory axial velocity fluctuation. 
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The complimentary solution of equation (16) is 
U c - C 2 J 0 (z) + C 3 Y 0 (z) 

where z = /-T Xr and X = ./— . 

v v 

The particular solution is obtained by the method of 
variation of parameters (Wylie [7]). 



U p = u l y l + u 2 y 2 
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The denominator of and u 2 ^ is the Wronskian of (y^,y 

(Abramowitz and Stegun [8]) 
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u 2 " = j q(z) R(r,x,t) 
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Entering the formulations for u and v into the expression for 
R(r,x,t) and then substituting the resultant expression 
equations (23) and (24) yields 



i(f> 
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(25) 
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where the various constants are listed in Table I 
Similarly 
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After integration to obtain the expressions for u^ and and 
substitution back into equation (18) , the complete expression 
for u is obtained. 



u = C 2 J 0 (z> + - X^* 
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The constant is made equal to zero by the boundary 
condition that u must be everywhere finite. The constant 
is solved for by the boundary condition that u = 0 at 
r = r Q . 
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The difficulty involved with this expression is obvious. 
Thus a second attempt at a theoretical solution was made to 
obtain a more manageable solution for the ease of calculation. 



B. SOLUTION # 2 OF PULSATING FLOW EQUATIONS 



In this second phase of the theoretical analysis, the 
assumption that v is independent of x is made. Thus, 
equation (4) becomes 
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1 8v 
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(28) 



Again, an oscillatory form for v (now independent of x) 
is assumed, namely v = R(r) e ^ Nt+< ^ . Substitution of this 
formulation into equation (28) yields a first order Bessel 
equation 
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(29) 



with the solution 



R(r) = (i^^ A r) (30) 

Incorporating this new expression for v alters the results 
for u^ and of the particular solution of u. Thje new 
expression for u^ becomes 
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As J v (x) and Y^ (x) are both cylinder functions, the 
expression for can similarly be expressed by substituting 
Y^ for J^. The constants X 2 j are tabulated in Table II. The: 
final expression for u thus becomes 



u = C 2 J Q (z) + u 1 Jq(z) + u 2 Yq ( z) 



(32) 



The immensity of u has been somewhat reduced. However, 
it was felt that additional assumptions would yield a very 
workable solution. 

C. SOLUTION # 3 OF PULSATING FLOW EQUATIONS 

After obtaining the two previous solutions for u, a more 
simplified expression was still desired. An order of magni- 
tude analysis for u near the wall and for u near the center- 
line was made, these being the two main areas of interest 
during the investigation. These order of magnitude determin- 
ations did lead to the finding of a simplified solution 
applicable across the entire cross section of the pipe. 

The expression for v is 

v = C 1 e i(Nt+<f,) J 1 (z) (33) 
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Based on the boundary condition that v = 0 at r 



R C 1 e ;i - (Nt+<J,) J 1 (i 3/2 \r Q ) = 0 
which is satisfied either by = 0 or by 

R e i(Nt+ « J l( i 3 / 2 Ar 0 ) = 0 
Neglecting the trivial solution, then 
R e i(Nt+ *> J l( i 3 / 2 Xr Q ) = 0 
which leads to 
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or simplified 
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(37) 



(38) 



where x^ = Ar^ . This form suggests a phase angle dependence 
upon both time and frequency. For t=0 , the dependence of 
phase angle upon frequency of oscillation is- shown in Figure 
1. The constant in this eigenvalue formulation is 
indeterminate because of the linearity of the problem. 

However, as in all such cases, it is the form rather than the 
amplitude of the perturbation velocity which is of interest. 
With <p thus determined from these eigenvalue considerations, 
v is of the form v = e^ (Nt+<f>) where the arbitrary 

constant has been set equal to unity for calculation purposes . 
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Determinations upon u arise from equation (7) 
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Near the wall, the velocities u and v appearing on the right 
side of equation (39) approach zero. The velocity gradients, 

and , though large, are finite. Thus, near the wall, 
equation (7) becomes that treated by Uchida, namely, 



V V 2 a - i |£ - |H = 0 

p 9x 9t 



(40) 



The resultant solution of (40) with x = Ar is 
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The solution of this non-homogeneous Bessel equation is 
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In this expression for u, the last term is negligible 
compared to the first two. Thus, a complete solution for u 
applicable to the entire pipe cross-section is 



u = 



P N Ve 



int 



iN 



J Q (i 3//2 x) 

J 0 <i 3 / 2 x 0 > 



~ 1 



(44) 



This last form is equivalent to Uchida's final result for 
the axial fluctuation velocity component. Additionally, the 
companion normalized (C^ = 1) radial fluctuation velocity 
component is 

~ = e i(nt+<£) J (i 3/2 x) (45) 
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Ill . 



EXPERIMENTAL ANALYSIS 



A. EXPERIMENTAL APPARATUS 

The experimental trials were conducted in a horizontal 
Lucite tube 1.5 inches I.D. and 45 feet (300 diameters) long. 

A centrifugal fan was used to pull air from the stagnation 
conditions in the plenum chamber through the pipe and 
pulsator assembly. The inlet chamber consisted of a rectang- 
ular box three feet square in cross-section and four feet in 
length; the inlet chamber was fitted internally with five 
6-24 mesh wire screens to eliminate inlet turbulence effects. 
Joining the plenum chamber and the tube was a two foot long 
settling chamber. The well rounded inlet horn was a quarter- 
section of a 9 inch by 18 inch ellipse molded of fiber glass 
and polished smooth. The tube, rectangular inlet chamber, 
and settling chamber were leveled and aligned to eliminate 
any curvature effect. 

The phlsator assembly consisted of a sinusoidal oscillating 
piston arrangement driven by a MB 1250 Exciter. The pulsator 
was located at the end of the Lucite pipe immediately prior 
to the centrifugal suction fan. The piston arrangement 
consisted of a 12 inch diameter disc reciprocating over an 
11 inch diameter opening to the atmosphere within the 13 inch 
cubical outlet chamber. The pulsating pressure signal was 
created by vibrating the piston in a direction perpendicular 
to the direction of the main flow velocity. 
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Velocity measurements were performed with a 1248-L-T 
1.5-6 Thermo-Systems Inc. type cross-wire multi-sensor probe 
and associated electronics. The hot wire probe could be 
inserted at measuring ports located at 16, 20, 32, 40, 48, 

56, and 272 diameters downstream of the entrance. 

In all cases, the hot wire probes were mounted in a 
transverse micrometer mechanism, which was capable of placing 
the sensor to within .001 inch of any desired radial position 
except at the wall. Probe configuration limited placement 
at the wall to within .050 inch from the rigid wall. 

B. EXPERIMENTAL PROCEDURE 

The developing mean velocity profile was monitored at 
existing ports in the Lucite pipe located at 20, 40, 56, and 
272 diameters downstream of the entrance. The quantitative 
data was recorded at these four stations for Reynolds Number 
values (R = U^^r/y) of 2848 and 5398 respectively. 

The fluctuating pressure differentials for high and 
medium mass flows were obtained using BRUEL & KJAER i inch 
condenser microphones (Model 4136) . The mean pressure 
differentials were obtained using micro manometers. 

To measure the pulsating velocity components, the cross 
wire probe was installed at station 1, twenty diameters 
downstream of the entrance. Radial distribution of velocity 
amplitudes and phase shifts were taken at sinusoidal pressure 
oscillations of 15 cps and 20 cps and at high and medium 
mass flow rates. All time-dependent data were recorded on a 
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two-channel Sanborn Model 296 recorder. The physical setup 
of the electronic equipment utilized to measure the oscillating 
velocity components and pressures is shown in Figure 2. 

On the two-channel recorder v and u data was recorded 
alternatively; the input pressure signal was continuously 
monitored . 

C. EXPERIMENTAL RESULTS 

Initially, mean flow velocity data was evaluated at station 
1 where the oscillating velocity determinations were to be 
obtained. A smooth curve was fitted to the mean flow data by 
the least squares method. Figures 3 and 4 illustrate the 
computed plots for flow Reynolds Numbers of 2848 and 5348 
respectively. Table III gives the numerical results for the 
developing velocity profile at the succeeding stations. 

Secondly, verification of the phase shift, <f> , inherent in 
the radial velocity expression was sought. After choosing 
an arbitrary t = 0 setting upon the records of the pressure 
and velocity signals, a subsequent marking of parallel times 
allowed for a direct verification of analytical values as 
shown by the representative tracings in Figures 13 and 14 
for a pulsation frequency of 20 cps . While equipment limita- 
tions forbid a conclusive verification at this frequency; 
nevertheless, the appearance of the 90 degree phase shift 
between the axial velocity and pressure signals on the center- 
line which is evidenced in Figure 15 and which is also 
predicted by theory makes the radial results less tentative. 
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Finally, a comparison of theoretical and experimental 
values for both radial and axial fluctuation velocity 
components was determined as shown in Figures 5-12. In the 
case of the axial velocity component, the comparison between 
experimentctlly measured values and theoretically calculated 
values is direct whereas for the radial velocity component 
the theoretical curves incorporated an arbitrary amplitude 
factor . 
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IV. CONCLUSION 



The final theoretical conclusions follow from the anamolous 
physical situation represented by equation (28) which demon- 
strates the viscous decay of the radial velocity component. 
Removing the x dependence of v, removed the production term 
— 3 v 

(u which was necessary to generate the v term initially. 

Thus, subsequent theoretical determinations are made upon a 
velocity component which is not strictly allowed by the linear 
analysis pursued. This anamoly is represented by the 
appearance of the arbitrary constant in the final expression 
for the radial velocity component. 

A second manifestation of this anamoly is evident in the 
independence of time that exists in the final expression for 
v since 

(nt+<}>) = tan ^ ber^ x^/bei^ x^ 

Thus the initial form of v which follows from this result must 
persist. Both the form and persistence (since RMS evaluations 
are involved) seem to be suggested by the comparison of the 
theoretical and experimental values for the radial velocity 
component portrayed in Figures 9-12. Again it is to be noted 
that the theoretical curves here incorporate an arbitrary 
amplitude factor. However, as in all linear eigenvalue 
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problems, it is the form rather than the absolute value of 
the eigenfunction which is important and which is verified 
experimentally here. 

The comparison of values in the case of the axial velocity 
component as expressed in Figures 5-8 is more conclusive. 

This data, though limited in nature verifies previous (Richard- 
son [9]), theoretical and experimental determinations and 
thus helps substantiate the parallel radial data which was 
simultaneously recorded. 

These mean determinations portray another characteristic 
of pulsatory flows, namely, that the maximum mean flow does 
not occur at the center of the pipe but in a region near the 
wall. The location of this region approaches the wall as 
pulsation frequency increases . As evidenced in both the 
radial and axial data, the high value of the argument in the 
Bessel functions involved in the velocity expressions makes 
the wall area values predominant. Thus the central portion 
of the flow is recognized as being practically unaffected by 
pulsation. Moreover, the value of axial fluctuation velocity 
is a non-zero constant on the pipe axis - another result 
verified by the data in Figure 5-8. In contrast, the zero 
value of radial velocity fluctuation on the pipe axis is 
evidenced by Figure 16 which portrays the radial velocity 
component recorded with the same gain readings as pertain to 
all other tracings shown. The appearance of this figure is 
that of a filtered noise signal and demonstrates the favorable 
signal to noise ratio evidenced by the other recorded data. 
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In addition to the results concerning the existence of 
a radial fluctuation velocity component, the second new 
result warranted by the accompanying analysis was the appear- 
ance of a frequency dependent phase shift in the formulation 
of this velocity component. Results concerning this phase 
shift are shown in Figures 13 and 14 where determinations 
for t = 0 values approach those values of phase shift 
suggested by the analysis for a pulsation frequency for 20 cps. 
Equipment limitations did not allow verification of non-zero 
time values. However, as stated previously, the appearance 
of a 90 degree phase shift in the axial velocity and pressure 
data recorded in Figure 15 tends to substantiate the companion 
radial velocity data. 

In conclusion, the limited experimental program reported 
here only demonstrates trends toward the parallel theoretical 
determinations. A more comprehensive experimental program 
would require both new electronic filters with band pass 
capabilities in the .01 cps range and a new pressure monitoring 
device to replace the microphone system which has severely 
attenuated response below 20 cps. In view of these require- 
ments, a re-evaluation of analytical results is recommended 
before a more extensive experimental phase is commenced. 
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TABLE I 



Constants for series solution of v; solution # 1 of the 
Pulsating Flow equations. 
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9 * X 9 = ( - 5b 0 B 8- 4b l B 7- 3b 2 B 6 ) 

V A 
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TABLE II 



Constants for v independent of x solution; solution # 2 
of the Pulsating Flow equations. 



1. 



'21 



c .u n iT 

4 0 

r„ ^ /nv 



2 . 



X 22 



C « V ^ 

4v 




1/2 



3. 




4r 0 2 /nv 



4. 



'24 



C 4 U 0 1T>/ri 
2 v 




5. 




c 4 y/=r 

2Tq^ /irv 
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TABLE III 



Coefficients of the polynomial approximations to the velocity 
profiles of the form: 



u 

U M 




2 



+ b. 




3 



x/d 


mass flow rate 


% 




^2 


^3 


20 


high 


0.9915 


-.0964 


1 . 3958 


-2.4025 




medium 


0.9760 


.3159 


-.1252 


-1.2745 


40 


high 


1.0137 


-.4726 


2.3368 


-2.9733 




medium 


0.9804 


-.1590 


0.5294 


-1.7836 


56 


high 


1.0010 


-.1717 


1.3127 


-2.1629 




medium 


0.9996 


-.1579 


0.9104 


-1.8154 


272 


high 


0.9653 


.6795 


1.6906 


-0.0106 




medium 


0.9904 


.1320 


-1.5630 


0.4163 
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Phase Shift (degrees) vs Frequency (cps) of 
v compared with pressure signal at t = 0. 
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Figure 2. Physical arrangement of electronic 

equipment for measuring u, v, and p. 
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Figure 3. Developing Velocity Profile; x/d=20, 
R=2848. 
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Figure 4. Developing Velocity Profile; x/d=20, 
R=534 8 . 
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Figure 5 . 



r/r o 

Axial Velocity Component; R=5348, 
Freq = 20 cps . 
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Figure 6. Axial Velocity Component; R=5348, 
Freq = 15 cps 
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Figure 7 . Axial Velocity Component; R=2848, 
Freq = 15 cps . 
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Figure 8. Axial Velocity Component; R=2848, 

Free = 15 cps . 
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Figure 9. Radial Velocity Component; R=5348, 
Freq = 20 cps . 
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r/r Q 

Figure 10. Radial Velocity Component; R=5348, 
Freq = 15 cps . 
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Figure 11. Radial Velocity Component; R=2848, 
Freq = 20 cps . 



42 




Figure 12. Radial Velocity Component; R=2848, 
Freq = 15 cps . 
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Figure 13. Recorder Tracing of p (left) and 
v (right) Radius = 0.7, R=2848, Freq = 

20 cps , <(> ~ 45° . 
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Figure 14. 
v (right) 
20 cps , 4> 



Recorder Tracing of p (left) 
Radius = 0.7, R=5348, Freq = 
~ 45°. 



and 
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Figure 15 . 
u (right) 
20 cps , 4> 



Recorder Tracing of p (left) and 
Radius = 0.0, R=5348, Freq = 

~ 90° . 
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Figure 16. Recorder Tracing p (left) and 
v (right) Radius = 0.0, R=5348, Freq = 
15 cps , v * 0 . 
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